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Abstract 

We discuss the behavior of response functions in phase ordering kinetics 
within the perturbation theory approach developed earher. At zeroth order 
the results agree with previous gaussian theory calculations. At second order 
the nonequilibrium exponents A and \r are changed but remain equal. 
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I. INTRODUCTION 



The scaling properties of response functions for phase ordering systems has been the sub- 
ject of some recent [1] interest. We study these properties here using the ^-expansion method 
developed previously [2] to extend perturbation theory beyond the gaussian level. We find, 
while there are no corrections generated at second order for the exponent a governing the 
response function, there are second order corrections for the nonequilibrium exponents A and 
Xr and X — Xr. We also find the associated scaling function at second order in perturbation 
theory. 

If ip{r,t) is the scalar ordering field whose dynamics are driven by a time dependent 
Ginzburg Landau (TDGL) model (described in detail below), then we can study the corre- 
lation function 



This quantity has been studied in a growth kinetics context [3] for a variety of systems 
and shows the scaling behavior C{r,ti,t2) — F{r/L{ti),ti/t2) where L{ti) is the charac- 
teristic growth law L{ti) ~ t^^^ and the growth exponent is ^ = 2 for the nonconserved 
order parameters studied in this paper. Focus has been on the onsite correlation function 
C{0,ti,t2) = F{0,ti/t2). For ti > t2 we have F{0,x) ^ x'^l^ where A is the well studied 
[4,5] nonequilbrium index. 

We can also introduce an external field, i^(l), conjugate to the order parameter and 
define the response function 



evaluated at zero field. For the local (ri = r2), response function it has become customary 
to write, > ^2, 



C(ri - r2, ti, t2) = (V'lri, ii)'0(r2, ^2)) 



(1) 




(2) 



X(0,ii,i2) 



(l+a) 



/(^lA2) 



(3) 



where for large x — ti/t 



'2, 



2 



f{x) ^ . (4) 

The goal in the analysis is to find the exponents a and and the scaling function f{x). 
Here we focus on the regime where both ti and t2 are large enough for the system to be in 
the scaling regime. 

Let us review some of the results found previously. The exact result [6] [7] for the one- 
dimensional Ising model gives a — 0, f{x) ~ (x - and Xr = 1. The solution of the 

problem in the large n limit [8] gives the exponent 

a = (d-2)/2 , (5) 

the scaling function 

and Xji = d/2. Berthier, Barrat and Kurchan [9] have carried out a gaussian auxiliary 
function approximation calculation for this problem. This has been extended to treat one 
dimension in [10] with the results: 

a = (d-l)/2 (7) 

and f{x) is again given by Eq.(6). Henkel, Pleimling, Godreche and Luck in [11], using 
conformal invariance methods, have derived a form of the scaling function: 

This is the same form as Eq.(6) if Xr/z = d/A. 

Numerical work [12] has focussed on smoothed integrals of the fundamental response 
function ^2)- This smoothing procedure helps with numerical sampling but can lead to 
qualitative differences between the fundamental response function and its smoothed counter 
parts for lower dimensions. 

We focus here on the computation of the local response ^2) = x(0, ^2) to second 
order in the 0-expansion developed in ref. [2] for a scalar order parameter. At zeroth order 
we find 



Xo(ii,i2) =t2"'~"/o(^iA2) , (9) 
with the exponent a given by Eq.(7) and 



/ 2 x'^l'^ 

where the constant is defined in section 4 below. 

At second order in perturbation theory we again find the scahng form given by Eq.(3) 
with a = \{d—\) but the scahng function is given now by 

/ 9 7.'»'-i/2 

where uj is determined by Eq.(14), v is given by Eq.(166), Arf(a;) is given by Eqs.(140) and 
(136) and g{x) is given by Eq.(157) and (163). In the large x limit, where lS.^{x) and g{x) 
go to a constant, we can use Eq.(ll) to identify the nonequilibrium exponent 

\R^d^\-2u-2v . (12) 

Let us use the same self consistent procedures for evaluating Ai? as used in I for A and 
a;. We found in I, at this same order, that 

A=2+-^337ST (13) 

with CO given as the solution to 

2a, + a,V(/Q+^) = l + ^ , (14) 

Jo [1 + zY 2 T{d) ^ ' 

and Kd is given by Eq.(164). Solving for to from Eq.(14) and inserting this into Eq.(13) gives 
the values of A shown in table 1. Notice, however, that if we eliminate Kd from Eq.(14) using 
Eq.(166) for v and eliminate in Eq.(14) in favor of X — d/2 using Eq.(13), we obtain 

2uj + 2v + {X- d/2) = 1 + ci/2 . (16) 



This can be rewritten as 



X^l + d-2u;-2v 



(17) 



Comparing with Eq.(12) and we find 



A = A 



R 



(18) 



and, as was found at lowest order, the nonequihbrium exponents are equal. These cancella- 
tions, coming from two very independent calculations of A and Xr, serve as a severe check 
on the validity of the algebra carried out in each. 

In the next section we describe the setting up of the perturbation theory used to obtain 
these results. 
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TABLES 



dimension 


X = Xr 


u 


1 


0.6268.. 


0.4601.. 


2 


1.1051.. 


0.6877.. 


3 


1.5824.. 


0.9067.. 


large 


d/2 


i(V2-l) 



TABLE I. Second order values for the exponents A = Ar and the parameter uj from ref.( [2]). 
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II. INTRODUCTION OF ORDERING AUXILIARY FIELD 



Consider a system where the ordering kinetics are driven by the simplest nonhnear time 
dependent Ginzburg-Landau (TDGL) model. Assuming that we quench a system from a 
disordered high temperature state to zero temperature, where the thermal noise is set to 
zero, we have the equation of motion. 



dil) 



r(-V'{iP)+cV^iP + B)+5{ti-to)ipo ■ (19) 



dt 

where ^(■0) is the driving potential. Typically this is chosen for simphcity to be of the ip"^ 
type: V = |(1 — i/j'^)'^. We have also included an external field S(r,t), conjugate to the 
order parameter, in the equation of motion. Choosing units of time and space such that 
r = c = 1 we can write 

A(1)V'(1) = -V'ii/Jil)] + B{1) + 5(h - (20) 

where the diffusion operator 

A(l) = |^-V; (21) 

is introduced along with the short-hand notation that 1 denotes (ri, ti). 

We show, extending the analysis in I to include an external field, that the order parameter 
■0 can be divided into an ordering component a and an equilibrating component m, 

-0 = cr[m] + u[m] (22) 

where cr[m] is the solution to the Euler- Lagrange equation 

^ a. = V'[a[m\\ (23) 

connecting the degenerate states = ±00, ^'(±'0o) = 0. In this equation m is taken to be 
the coordinate. In the case of a 0'*^ potential we have the solution 

cr[m] = tanh{-j=m) . (24) 



It is shown in ref. [2] that the equihbrating field u decays exponentially to zero at long times. 
The scaling properties of the theory are carried by the ordering field a[m]. 

It was shown in I that the theory is self- consistent if the auxiliary field saisfies an equation 
of motion of the the form 



where S is a function of m which must scale as ~ i/L{t) in the scaling regime and self- 
consistently generates ordering. In particular S must be such that (m^(l)) ^ L'^{ti). In I 
and here we study the nonlinear model where 



and, if scaling is to hold, we find self-consistently that i{ti) ~ 1/L{ti) and is independent 
of the field m. 



Let us consider the field theory associated with the equation of motion for m(r, t) given 
by Eqs.(25) and (26). Our development will follow the standard Martin-Siggia-Rose [13] 
method in its functional integral form as developed by DeDominicis and Peliti [14]. In the 
MSR method the field theoretical development requires a doubling of operators to include 
the field M which is conjugate to m. We give here an overview of the development in I 
needed here to treat the response function including the coupling to an external field. 

Following standard procedures, described in more detail in I, averages of interest are given 
as functional integrals over the fields m, M and weighted by the probability distribution 
P[m,M]: 



A(l)m(l) = S(l) + B{1) + S{ti - io)mo(ri) 



(25) 



S(l)=e(ti)a(m(l)) 



(26) 



III. FIELD THEORY FOR AUXILIARY FIELD 




(27) 
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and 

Z{h, H)^ j VrnVMe^^^""'^^ . (28) 

The action takes the form 

Arim, M) = A{m, M) + J dl [h{l)m{l) + H{1)M{1)] (29) 



where 



A{m, M) = -iJ dlM{l) [A(l)m(l) - ^il)a{l) - B{1)] 

-^J dl j d2 M(l)no(12)M(2) (30) 



with 



Tl^{l2) = 8{h-h)8{h-t2)g{v^-V2) . (31) 

In these equations we use the notation, J dl — J dtid'^ri , and assume that the initial field 
mo(r) is gaussian and has a variance given by 

< mo(ri)mo(r2) >= g{ri - rs) . (32) 

We can generate correlation functions as functional derivatives in terms of sources h and H 
which couple to the conjugate fields. 

The fundamental equations of motion are given by the identities 

^0^^^^""'^^^'"° ^^^^ 

<^^T(m,M)>,= (34) 

where the subscript h indicates that the average includes the source fields h and H. Taking 
the derivative in Eq.(33) we obtain 

z[A(l)<m(l) >h -Qi{l)] 



= - y d2 no(12) < M(2) >h +H{1) - iB{l) 
where the nonhnearities are included in 



(35) 



gi(i) = e(i)<f^(i) >h 



(36) 



Eq.(34) gives 



■i A(l)(M(l)) + gi(l) =/i(l) 



where 



A(l)^^ + V? 



and 



Q,{l)^mMm{l))M{l)) 



(37) 



Clearly we can go on and generate equations for all of the cumulants by taking functional 
derivatives of Eqs.(35) and (36). Let us introduce the notation that GAi,A2,....,A„{^'2---'n) is 
the n*'* order cumulant for the set of fields {^i, A2, An} where field Ai has argument (1), 
field A2 has argument (2), et cetera. This notation is needed when we mix cumulants with 
m and M. As an example 




< m(4) >h 



(38) 



SH{l)Sh{2)Sh{3) 



As a short hand for cumulants involving only m fields we write 



< m{l) >h 



(39) 



5h{n)5h{n - I) ■ ■ ■ 5h{2) 



The equations governing the n order cumulants are given by 



I A(l)G'Mm...m(12...n) + g,(12...n) =0 



(40) 



and 
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[A(l)G'„(12...n) - g„(12...n)] ^ - J dl no(ll)G'Mm...m(i2...n) 



(41) 



The Q's are defined by 

Qn{l2...n) = 

Qn(12...n) 



n— 1 



Sh{n)6h{n - 1) ■ ■ ■ Sh{2) 



Qi(l) 



(42) 



STn-l 



5/i(n)5/i(n-l)-- -(5/^(2) 
With this notation the equations determining the two-point functions can be written as 



(43) 



A(1)G'm„.(12) + Q2(12) =5(12) 



i [A(1)G'2(12) - g2(12)] --Jdl no(li)G'Mm(12) . 



(44) 



(45) 



This is all formally exact. In order to develop perturbation theory the next step is to 
show that Qi{l) and Qi{l) can be expressed in terms of the singlet probability distribution 



Ph{x,l) =< S{x-m{l)) >h 



(46) 



One finds 



gi(l)= jdx ^{l)a{x) Ph{x,l) 



Qi(l) ^1 dx C{l)<^i{x) 



<M(1) >ft + 



5H{1) 



Ph{x,l). 



(47) 



(48) 



Then any perturbation theory expansion for Ph{x, 1) will lead immediately to an expansion 
for Qi{l) and Qi(l). We can then obtain Qn and Qn by functional differentiation. 

The development of a perturbation theory expansion for Ph{x, 1) begins by using the 
integral representation for the (5-function: 

/• dk 



/n k 



(49) 



where 



11 



(50) 



The average of the exponential is precisely of the form which can be rewritten in terms of 
cumulants: 



h, 1) = exp 



.8=1 



(51) 



where Gs(ll...l) is the s-order cumulant for the field m(l). 

Consider first the lowest-order contribution to Qn which does not vanish with the external 
fields h,H: 



Q,{12) ^ jdx e(l)<7(a;)^P,(x, 1) . 



(52) 



We have shown that the n order cumulants are of order | — 1 in an expansion parameter 
we will develop. Expanding $(A;, h, 1) in powers of the cumulants with n > 2 and keeping 
terms up to the 4-point cumulant, we obtain 



Ph{x, 1) 



^-^^3(111)^ + ^^.(1111)^ + - 



(53) 



where 



^r(^,i) 



dk 
2^ 



$o(/E,/i, l)e 



-ikx 



(54) 



and 



$o(A;,/i,l) = e^'=^i«e-^'='«^(ii) 



(55) 



Then, after taking the derivative with respect to h{2), setting the external fields to zero, 
and neglecting all cumulants with n > 2, we obtain 



^o{k,h = 0,1) = e-2 



(56) 



and 



g^°\l2) = Jdx ^{l)a{x) J ^e-'^''ikG2{12)e 



-ifc2S2(l) 
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where we have defined in zero external field 



S2{1) = ^2(11) =< m'(l) > . 



(57) 



In the scaling regime, where <S'2(1) is very large, we can replace a{x) — > ipo sgn{x) in the 
integral and obtain: 

(12) =e(l)V'oG'2(12) j dx sgn{x) J ^^fce-''=^e-5'='^^W 



e(l)V'oG'2(12) I dx sgn{x) 



In this case we integrate by parts in the integral over x and use -^sgn{x) — 25{x) to obtain 

gf (12) = e(l)V'oG2(12)2$o(0,l) 



= e(l)V^oG'2(12)i 



(58) 



7r52(l) 

Turning to (52(12) we note that it is given by taking the derivative of Eq.(48) with respect 
to /i(2). In the scaling regime where the characteristic length L{t) is large we can replace 
ai{x) — > '^q25{x). Then we need only consider 



Q2(12) = e(l)V'o / dx25{x) 



X 



G'm„.(12)+ < M(l) >h 



6 



6' 



Ph{x, 1) 



6h{2) 6h{2)6H{l)_ 

A key observation is that as we analyze contributions to Qn or Qn we will find that 
each term consists of products of correlation functions and response functions with legs tied 
together by factors defined by 

dk 



(j)p{l) = / dx sgn{x) I ^ik^P+'e-"'''^o{k, 1) 



27r 



(59) 



2 / ^ 
J 27r 



dS2{l) 



0o(i; 
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where we have used an integration by parts in going from the first to the second hne and 
defined 



,.(l)^2/|^e-i--«^^^ . (60) 

Each term in the perturbation theory expansion for Q„ or Qn will be proportional to factors 
of (pp. The perturbation expansion is ordered by the sum of the labels p on (pp. Thus a 
contribution with insertions 0i020i , each factor typically associated with different times, 
is of C(4). We refer to this expansion as the (p-expansion. It should be emphasized that at 
this stage that this is a formal expansion. At order n it is true that (pp a; L~(^^+^) which is 
small, however it will be multiplied, depending on the quantity expanded, by positive factors 
of L{t) such that each term in the expansion in (pp has the same overall leading power with 
respect to L{t). 

To see how this expansion works let us consider first the two-point quantity (52(12), 
defined by 

Q,il2) ^ J dx ai)i^oS9nix)^-^^ . (61) 

Using Eqs.(49) and (51) and taking the derivatives with respect to h{2), it was shown in I, 
in the case of zero external fields, that 

g2(12)=M(l)/ dxsgn{x)j |^e-*'=^$(A;, /i = 0, 1) 

xS(L^G..«(11...12) . (62) 
Since all odd cumulants vanish in the case of zero external fields we have 



00 (_lY]^2s 



(63) 



$(/c, /i = 0, 1) = exp 
where 

^2.(1) = G'2.(11...1) . (64) 
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Let us define the set of vertices 

Vp(l) ^ J dx sgn{x) j ^ik^P+^e-'^''^{k, /i = 0, 1) , (65) 

which can be written, after following the same set of steps in reducing the original expression 
for (f)p, to 

V,{l) = 2j'^k^^^Kh^Q,l) (66) 

which is independent of position. Then the quantity (52(12), which appears in the equation 
of motion for G2(12), is given in the form 

-i: 



g2(i2) = M(i) E 7^— V^v,(i)G2.+2(ii...i2) 



(67) 



where we have used the definition of Vs(l) given by Eq.(65) in the last step. 

It should be clear that the vertices Vs(l) are of at least 0{s) in the 0-expansion. By 
direct expansion of ^{k, /i = 0, 1) about ^^{k, /i = 0, 1) we obtain 

V.(l)=0.(l) + ^</>.+2(l)-^0.+3(l) + --- . (68) 

It was found self-consistently in I that order cumulants, hke ^'^(l), are of 0(| — 1). 
The terms in the expansion for Vs, given by Eq.(68), are of C(s), 0{s + 3), and 0{s + 5) 
respectively. 

Let us turn next to (^2(12). It was shown in I, in the same limit of zero applied field, 
that the nonlinear contribution to the equation of motion for Gmu can be written as 

oo / -i\s 

Q,{12) = e(l)V^o E ^V.(1)G(J::.1m.„(11...112) . (69) 

Before going on to discuss the perturbation theory calculation of the physical response 
function let us made sure the theory is sensible at zeroth order where, from Eqs.(67) and 
(69) 

gf (12) = e(l)V'o0o(l)G'2(12) = 6^0(1)^2(12) (70) 
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(12) = e(l)V'o0o(l)G'Mm(12) = u;o{l)GMmm 



(71) 



where for a scaling solution 



Uoil) = e(l)V'O0o(l) 



(72) 



must fall off as 1/ti for large ti. 



IV. GAUSSIAN THEORY 

Inserting Eqs.(70) and (71) into Eqs.(44) and (45) we obtain the equations for the re- 
sponse function: 



and the correlation function 



A(l)+u;o(l) g£L(12) = 5(12) 



(73) 



z[A(l) -(^o(1)]G(°)(12) = -J dlno(li)Gff„(l2) 



(74) 



It is not difficult to show that 



(75) 



Using this last result we have that the correlation function can be written as 

Gf\l2) ^jdlj d2iG(;L(ll)^Gi(22)no(12) (76) 

where no(12) is given by Eq.(31). 

The first step in the construction of the solution to these equations is to Fourier transform 
Eq.(75) over space: 

d 



+ q -uJo{ti] 



iG%{q,ht2) = 6{ti-t2] 



(77) 



This first-order differential equation has the solution 



iG^mMi<l,ht2) = 0{ti - t2)exp 
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^ ' dr (-q^ + a;o(r)) 



(78) 



and we have defined 

R{t,,h) - e^'^^ "^^'"-^^ . (79) 
Taking the inverse Fourier transform we obtain 



Let us turn our attention to the correlation function: Taking the Fourier transform and 
inserting the results for the propagators and Ho we obtain 

G^2\(lMt2) = e{h - to)e{t2 - to)R{h,to)R{t2,to)e-^''"^g{q) (81) 

where g{q) is the Fourier transform of the initial correlation function and T = — to- 
While we are primarily interested in the long-time scaling properties of our system, we can 
retain some control over the influence of initial conditions and still be able to carry out the 
analysis analytically if we introduce the initial condition 

m = goe-"^^"'^" (82) 



or 

Inserting this form into Eq.(81) and doing the wavenumber integration we obtain 

Let us turn now to the quantity R{ti, ^2) defined by Eq.(79). We assume that ujq has the 
form given by 

^o(l) = r^ ' (85) 

tl + tc 
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where tc is a parameter (or function of t which goes to a value tc for large times) such that 
the correlations of m has a smooth early time behavior, a; is a constant we will determine. 
Evaluating the integral 



I dr (jJo{t) = I dr — to In (———-\ , (86) 

Jt2 Jt2 tc-\-T \t2 + tc/ 



Rit^M = {P-^] ■ (87) 



we obtain 



.t2 + tc 

Inserting this result back into Eq.(84) leads to the expression for the correlation function 

Gf(r,M.)=,(o)r^r7^r;:^^ . (88) 

If we are to have a self- consistent scaling equation then the autocorrelation function (r = 0), 
at large equal times ti — t2 — t, given by 

:'(0),'+^ _ ^2oj-d/2 1 9o 



S'^\t)^e^-''"- , , (89) 



must have the form S2^\t) — AqI. Comparing we see the exponent ou must be given by 

^ = ^(1 + ^) (90) 

and the amplitude by 

(to+tcy- (8^)V ■ (^1) 

The general expression for the correlation function can be rewritten in the convenient 
form 

(r, hh) = V4°)(t0^f(t2)$o(tit2)e-^^/(^^+^^) (92) 

where 

d/2 



T + t^ + to 
18 



The nonequilibrium exponent is defined in the long-time hmit by 

A 

(94) 



^^(0,^2) _ (^{U+Qih+tc^ 



po{h)So{h) \ T + to + t, 
and we obtain the OJK [15,3] result 



■ (95) 

Looking at equal times we have that 

/.(-) = = (96) 

where the scaled length is defined hy x — r/At . fo{x) is just the well known OJK result for 
the scaled auxiliary correlation function. The connection to the physical order parameter 
correlation function is discussed in I. We will not need these results here. 



V. NONLINEAR RESPONSE 
A. General Expansion 

We are interested in the nonlinear response function 

= i((7(l)M(2))^,B=o 

= ^^(-(l))..^o (97) 

where the second line follows from Eqs.(27) and (30). We have now set B{1) — and we 
can express the reponse function in terms of a functional derivative of the singlet probability 
distribution function just as in the computation of Q2- 

r 
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The perturbation theory expansion for Ph{x, 1) was discussed in section 3. Using Eqs.(49) 
and (51) we have 

dx a{x) j — e-'^^$(fc, h, 1) J2 ^a,M(11...12) . (99) 



We can then set the source fields to zero and obtain 



-1 



(100) 



where we introduce the same vertices as in section 3. We have the exhcit expansion for 
the physical response function to the lowest two orders: 

X(l, 2) = z0o(1)G'^m(12) + ^0l(l) tllG'„^(1112) + . . . . (101) 

The leading term says that the physical response contains a term proportional to the MSR 
response function G^m(12). 

B. Zeroth Order 

At zeroth order we have the response function 

Xo(1,2)=z0o(1)gSL(12) . (102) 

0o(l) is given by Eq.(60) and iG^j^{12) is given by Eq.(80). Putting these together we 
obtain 



If we focus on the on- site response function: 

xo(h,t2)^xo(o,ti,h) ^ ^=J==e(h-t2) (^^)" .... ^ . ■ (104) 
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If we assume ti > t2 and write S2^\l) — Aot, we have 



/ TlAoti ' 






\t2 + tc) 



It is conventional to write this in the scahng form 

Xo(h,t2)^t^'-''Mx) (106) 
where x — ti/t2. Clearly we can identify l + a = l/2 + ci/2, or 

a = ^ (107) 

and the scaling function is given by 



and, using Eq.(90) for lu, we find in the large x limit the form given by Eq.(4), with 

Xr = d+l-2uj = d/2 = X . (109) 

C. Second Order 

At second-order we have two contributions. The second-order contibution to GmM can 
be read off from Eq.(165) in I with the result 

xi'^(12)=i0i(l)G'i^L(12) (110) 

where 

g21(12) = / dld2G%{ri)J:^'\l2)G%{22) (111) 
and the associated self-energy is given by 

E(2)(12) = ^ (-ia;i(l)) G%{12) (cf (12))' {-iM2)) . (112) 

cui is given by 
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Ml) 

The second contribution at second order is given by 



-1(1) = ^ • (113) 



X?'(12) = i0,(l)i^G<2LM(1112) . (114) 

We therefore need the lowest order expression for G^^^^j^ (1112) . The lowest order expres- 
sions for all of the various four-point cumulants were worked out in 1. g|||]„^j\^(1112) is given 
by Eq.(155) in I: 

gSLm(1112) = /cii3GSL(ll)(G^°^(ll))'(-^^i(l)Gi(i2) . (115) 
This second-order contribution is given then by 

x?\i2) = f ciiV^o|£:^|iG'(^L(ii)(Gr(ii))V(i)^Gi(i2) . (iie) 

In order to carry out the integrals contributing to x*^^-* it is convenient to write 

2^] ^^^^^ 



where we introduce the notation 



and 

(i/2 



= ;T7^^ (118) 



(12) = R{t„to)R{h,to)9o (1^) e-C^y^^ (119) 



where 

1 



(120) 



2{U + tj) 

Let us first evaluate X2^^(12) which is somewhat simpler. Inserting these forms for the 
correlation functions and response functions into the integral for X2^'*(l, 2) we obtain: 



(1, 2) = ^o^^^(ti - t2)R{h, t,)R\t,,to)glK{l2) (121) 
(bo 
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where 

K(U) = £ (I^f (l^f AU) (122) 

and the spatial integral is given by 

J(12) = / d'^fie V 2 J lie V 2 ^ 12 . (123) 

The spatial integral can be evaluated by completing the square in the gaussian with the 
result: 



j(12) = (_) e V ^24) 

where 

ao = ^'ii + ^'12 + 2aii . (125) 
Let us focus on the onsite correlation function (r^2 — 0) 

.(M.) ^ £ *.^(., (l;!)^'^ (^g {^^f (.iFTi^)" 

(1 \ 3ci/2 
— j dtii?'(ti,to)^i(i)<^' (126) 



where 



^il + bi2 + 2aii 

= [8(ii + fi) ((ii - t2){h + ii) + 2(fi - t2)(Mi))]"' ■ (127) 
We have from the zeroth-order theory 

sf\l)=Aoh (128) 

and 



^i(l) = 772 (129) 
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where u is given at zeroth order by Eq.(90). We have then 

1 /I \3'^/2 



UJ 1 / 1 \ ' 



where 



W = 8W0 (132) 



Using the result 



go^A^t^i^ny/^ (133) 



which follows from Eq.(91), and collecting coefficients, we have 

X?^ (ti, ^2) = - ^2)^(^1' i^)i"iK{t,, h) . (134) 

Looking at the time integral we change variables from ttot — yti and introduce s — t2/t\ 
so that 



^(^1, ^2) = f/^ [y^A m - «)(1 + y) + 2tl{l -y){y- s)] 



y V + y, 

= q'^n{s) (135) 

where 

<s) =1^^ [iT^y -s){l + y) + 2(1 - y){y - s)]-'/' . (136) 
Then the second-order correction is given by 

= ^000(1)2^^(^1 - t2)Rih, t2)h'"\{s) . (137) 
The zeroth order contribution is given by 

24 



Xo{tu t2) = V'o0o(l)^(ii - t2)R{h, h) (47r(ii - ^2))"''/' ■ (138) 

(2) 

We can then combine the zeroth and X2 contributions to find 

Xiti,t2) = Xoih, t2) + X2\ti, t2) = XoihM) (1 + Ad(s)) (139) 

where 

Ad(s) = u2^-\l - sYI^k{s) , (140) 

and k{s) is given by Eq.(136). Analytically one can show Ad(0) ^ ^ as 0? — > 00, and 
Ad{s) ^ ^{1 - s) as s ^ 1. 

(2) 

We turn next to the evaluation of G Jj^. Inserting the response and correlation functions 
and using the properties of the R{ti, ^2) we have 

G^L(12) = - t^)R{t^, t,) £ dh £ df2f ^i(i)^i(2) 

X {bnbi2b-22f^ai-^R''{iuto)R^{h:to)JrnM (141) 

where 

/(jd-f^dd-f^ _6l 2 _b2 2 2 e~2'''ll 

where 

61 = &1I (143) 
62 = 612 + 2ai2 (144) 
h = &22 (145) 
V = 6162 + &2&3 + hhi (146) 

and 

25 



a ^ (b^' + b^' + b^'Y' . (147) 
Again restricting ourselves to the onsite response function 

gSL(0, h, h) = -^J^^(^iti - h)R{h, h) £ dh rft2|u;i(i)c^i(2) (148) 

><{hiibi2h-22f''42R\iiM)R\i2M^ ■ (149) 
After grouping the various multiphcative terms, using Eq.(133), gives 

cJ^mIO, tu t2) = -te{h - t2)R{h, t2)u;^j^^J{h, t2) (150) 

where 

Ah, t2) = r dh r dUiih?^-^ {wf^ (151) 

Jt2 Jt2 



with 



W= , ^''f' , , , . (152) 

1 + 62&3 + ^'2&r 



After some algebra we find 



W ^ - — ^ 1- ^— . (153) 

8{ti + t2)[{ti-t2){3ti-t2)-2{ti-t2)^] ^ ' 

After inserting this result for W back into Eq.(151), using a; = d/2 + 1, which is valid at 
this order, and letting ii — {ti — ^2)^/1 and ^2 = (^i — ^2)^/2 we obtain 

where 

dyi ry^dy2f ym \ '^''^ 1 

Jr Vi Jr Vo Ul+y2/ [3^1 - ^2 - 2(^1 - yo)"!''/" ^ ' 



9{r) 
and 



r = ^^ = — i . (156) 

tl-t2 h/t2-l ^ ' 

26 



A more useful form for g{T) is obtained if we let 1/2 — yiz to obtain 

™1+T dm /•! 



9(r)^ r dzf(z,y^) (157) 



where 



^^"'^^^^[(l + .)((3-.)-2|/,(l-.)2)]^/^ ■ ^^^^^ 
This means that the MSR on site response function can be written up to second order 

GmM{Q, h, h) ^ -ie{h - h) (1 + u'2'-'g{T)) . (159) 

Pulling together all of the results for the physical response function function we have up 
to second order 

X{h,t2)^t^'-''fih/t2) , (160) 



with 



a=l(d-l) (161) 



and 



w-1/2 



(162) 



where x — ti/t2- It is not difficult to show that in the small r limit 

g{T) = -KalnT + g{T) (163) 

where 

i^d also appeared in the analysis in 1. If our perturbation theory results are to make sense 
we must exponentiate the —K^ In r contribution: 

1 - t; T fa T"" = (x - 1)^ (165) 
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where in the last teerm we have used Eq.(156) and 

V = uj'^2'^-^Ka . (166) 

Putting this result back into Eq.(162) we obtain our final result for the scaling function 
given by Eq.(ll). 

VI. CONCLUSIONS 

We have used perturbation theory to explore the nature of the scaling solutions for the 
local scalar order parameter response function. At lowest, gaussian level, order we find 
agreement with previous work. Going to the next order we find that the scaling index a, 
defined by Eq.(3), is unchanged from its gaussian level value of a = (l/2)(d — 1). The 
non equilibrium exponents A and Xr are found to be equal at lowest and second order in 
perturbation theory with explicit values given as a function of d. Expressions for the scaling 
function are also available. 

It is not at all clear if there is a general proof that \ — \r. The validity of this result 
at second order in perturbation theory is very suggestive. It will also be interesting to see 
whether, within this perturbation theory, this result holds for the case [16] of vector order 
parameters (n > 1) where the analysis is more involved. 
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